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ABSTRACT 

O. 

^ . Because H 2 formation on dust grain surfaces completely dominates gas-phase 

H2 formation in local molecular clouds, it is often assumed that gas-phase for- 
mation is never important. In fact, it is the dominant mechanism in a number 
of cases. In this paper, I briefly summarize the chemistry of gas-phase H 2 forma- 
tion, and show that it dominates for dust-to-gas ratios less than a critical value 
T> cr . I also show that T> CI is simple to calculate for any given astrophysical situ- 
ation, and illustrate this with a number of examples, ranging from H2 formation 
O ! in warm atomic gas in the Milky Way to the formation of protogalaxies at high 

"c^ redshift. 
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1. Introduction 

In local molecular clouds, molecular hydrogen (H 2 ) forms primarily on the surface of 
dust grains: two hydrogen atoms are adsorbed onto the surface of the grain and react to 
form H2, which subsequently escapes back into the interstellar medium. However, H2 can 
also form in the gas-phase, primarily through the reactions 

H + e -» H-+7 (1) 
H- + H -» H 2 + e, (2) 

although some also forms via the slower reactions 



(3) 
(4) 
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In dense gas, three-body reactions can also be important (Palla, Salpeter & Stahler 1983), 
but these are ineffective at number densities n < 10 8 cm' 3 . 

Discussions of H 2 formation have tended to concentrate on the role played by dust, with 
little attention given to the gas-phase reactions. However, as I show in section 4, in some 
circumstances these reactions can dominate the H 2 formation rate. 

In this paper, I briefly outline the chemistry of H 2 formation and show that it is easy to 
identify a critical dust-to-gas ratio V cr , above which grain- catalyzed formation dominates. I 
illustrate the method by applying it to various situations of astrophysical interest, and show 
that it can be a useful tool for estimating the importance of gas-phase H 2 formation. 



2. The formation of molecular hydrogen 

2.1. Gas-phase formation 

Most of the molecular hydrogen that forms in the gas-phase does so via the formation of 
an intermediate H~ ion, as outlined in reactions 1 and 2 above. The first of these reactions 
occurs much more slowly than the second, and so the equilibrium abundance of H~ is small 
and is rapidly reached. Thereafter, the H 2 formation rate is determined by two factors: 
the rate at which H~ forms, and the fraction of H~ ions that survive to form H 2 . The 
latter quantity is determined by competition between H 2 formation via reaction 2 and H~ 
destruction by mutual neutralization with H + ions 

H + H+ -»• 2H, (5) 

and by photodetachment by the incident radiation field 

rT+7^H + e. (6) 

Various other reactions also destroy H~, but these are either significantly slower than those 
above, or become important only at high temperatures, in which case any H 2 that does form 
will very quickly be collisionally dissociated. For more details, the reader is referred to the 
recent reviews of Abel et al. (1997), Galli & Palla (1998), Stancil et al. (1998) and Lepp 
et al. (2002). 

If we assume, for simplicity, that H~ has already reached its equilibrium abundance, 
then we can write the H 2 formation rate as 



(7) 
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where Hi is the number density of species i, and where the rate coefficients hi for the various 
reactions are listed in table I. 

If H 2 formation via reaction 2 occurs much faster than the destruction of H~ by the 
other reactions, then this reduces to 

#h 2 ,h- - hncUH, (8) 
in other words, the H 2 formation rate is approximately the same as the H~ formation rate. 

If, on the other hand, mutual neutralization dominates over H 2 formation or photode- 
tachment as a means of removing H~, then equation 7 becomes 

Rh 2 ,h- - hn e n H ^, (9) 

where x = «h+ / n u is the fractional ionization of hydrogen. As long as n c ~ nn+ , this 
equation can be further simplified to 

#H 2 ,H- ^ -j— %• (10) 

Comparing this equation with equation 8, we see that for a small fractional ionization 
i?H 2 x i but that once the fractional ionization becomes large enough that mutual neutral- 
ization dominates, i?H 2 becomes independent of the ionization: although increases in x still 
increase the H~ formation rate, this is balanced by the increase in the mutual neutralization 
rate and consequent decrease in the fraction of H~ ions surviving to form H 2 . This change 
in behaviour occurs for fractional ionizations near a critical value x cr , defined by 

x CT = (11) 

The precise value of x cr is somewhat uncertain, due to the significant uncertainty that remains 
in the determination of the mutual neutralization rate. In this paper, I have chosen to adopt 
the rate listed in Galli & Palla (1998), which is derived from the data of Moseley et al. (1970). 
This is a conservative choice, in that it gives the lowest value of x cr ; other possibilities include 
the rates of Duley & Williams (1984), Dalgarno & Lepp (1987) and Croft et al. (1999), with 
the last-named being preferred by the most recent compilation (Lepp et al. 2002). For the 
temperature range of interest, the Galli & Palla rate gives us a value x cr ~ 5 x 10 -3 , with 
only a slight dependence on temperature. The alternative rates typically give values of x cr 
that are factors of a few larger. 

When the destruction rate of H~ ions is dominated by photodetachment, we obtain 
another limiting case of equation 7 

#h 2 ,h- ^ hneTiH-^- (12) 
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This can be written as the H formation rate divided by a suppression factor f ra( i- 

^h 2 ,h- ^ —7 , (13) 

/rad 

where 

/rad = T~~~ ■ (14) 

k 2 n H 

Determination of the photodetachment rate, and hence / ra( j, requires knowledge of the 
incident radiation field. Provided that the opacity of the gas is low, we can write the 
photodetachment rate as 

h = An ^±dv, (15) 
Ah hv 

where J u is the mean specific intensity, a u is the photodetachment cross-section (de Jong 
1972) 

a v = 7.928 x 10 5 ^ - cm" 2 , (16) 

and where hu th = 0.755 eV is the energy threshold for H photodetachment. 

Evaluating equation 15 for the local interstellar radiation field, as estimated by Mathis 
et al. (1983), gives 

k 6 = 4.2 x Hr 7 s-\ (17) 

and so, locally, 

/ rad = 3.2 x lO 2 ^ 1 . (18) 

If the opacity of the gas is high, then this will overestimate the effects of radiation. 
However, at frequencies near the H~ photodetachment threshold, the continuum opacity of 
interstellar gas is low and absorption is dominated by dust. Consequently, a high opacity 
implies a high dust content, in which case grain- catalyzed formation will dominate. 

These limiting cases provide useful insight into the physics of gas-phase H 2 formation, 
but in general we must use the full form of equation 7, which we can rewrite as 

Ru2 H ^ = klUcU}l (19) 

1 ~\~ X j X CT + / ra d 

If we now turn to H 2 formation via the ion, we find that the basic chemistry is 
remarkably similar. is created by the radiative association of H and H + (reaction 3), and 
destroyed by H 2 formation (reaction 4), dissociative recombination 

H+ + e -> 2H, (20) 
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and photodissociation 

H+ + 7 ^H + H+. (21) 

As with H _ , the formation of the molecular ion is the limiting step, with subsequent 
reactions occurring orders of magnitude faster. If we again assume that the abundance 
has reached equilibrium, then we can write the H 2 formation rate as 

Ru 2 ,m = k ^ n » hnH + Z e + h > (22) 

which has the same form as equation 7. Indeed, we can rewrite it as 

k 3 n R+ n u 

only now 

x CT = (24) 

and 

/rad = T~~~) (25) 

where the photodissociation rate, k 8 , is calculated in a similar fashion to the photodetach- 
ment rate above. Evaluating these, we find that x CT for is typically an order of magnitude 
larger than for H _ , and that in the local ISM 

/rad ^ ^H 1 , (26) 

where I have again used the Mathis et al. (1983) radiation field, together with the H^" 
photodissociation cross-section from Shapiro & Kang (1987). 

Comparing these values with those for H~, we see that is significantly more ro- 
bust. However, it forms at a much slower rate (between two and three orders of magnitude, 
depending on temperature) and so in most cases H~ dominates. Nevertheless, there are 
exceptions, as we will see in section 4. 

Finally, a few other possible mechanisms have been suggested for gas-phase H2 for- 
mation. Latter & Black (1991) propose that H 2 can form as a result of direct radiative 
association 

H + H -> H 2 + 7, (27) 

provided that one of the hydrogen atoms is in an excited electronic state. Rawlings et al. 
(1993) show that a more efficient mechanism is formation of by associative ionization 

H + H -> H+ + e (28) 
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with the H2 thereafter forming H 2 by reaction 4 above. This mechanism again requires one of 
the hydrogen atoms to be in an excited atomic state. However, this requirement means that 
in general these reactions are not important, as the necessary population of excited atomic 
hydrogen is only found in a few unusual circumstances, such as in the universe immediately 
after recombination. 



2.2. Grain-catalyzed H 2 formation 

Despite its importance in local interstellar chemistry, the rate of H 2 formation on dust 
grains is still uncertain. In local molecular clouds, observations suggest a formation rate 
(Jura 1975) 

kg ~ 3 x l(r 17 cm 3 s" 1 . (29) 

Observations of H 2 in the LMC and SMC with the FUSE satellite (Tumlinson et al. 2002) 
suggest a value that is an order of magnitude smaller, but this is consistent with the under- 
lying rate per unit dust mass being the same, since the mean dust-to-gas ratio within these 
galaxies is significantly smaller than in the Milky Way (Issa et al. 1990). 

Unfortunately, direct measurements of this kind can only give us information about H 2 
formation in physical conditions that are easily accessible to observations, and provide little 
basis on which to predict the H 2 formation rate in different regimes. For this, we must turn 
to theory. 

A large body of theoretical work exists on the subject of H 2 formation on grains (see, for 
example, the review of Pirronello et al. 2000, and the many references therein), stretching 
back almost forty years to the pioneering work of Gould and Salpeter (1963). In a highly 
influential paper, Hollenbach et al. (1971) parameterized the H 2 formation rate as 

R = 0.5v H a d Sf a n K n d , (30) 

where n d is the number density of dust grains, a d is their mean geometric cross-section, % is 
the mean velocity of the hydrogen atoms striking the grains, S is the sticking coefficient (the 
probability that a hydrogen atom striking the grain will stick to the surface) and f a is the 
fraction of adsorbed hydrogen atoms that actually form H 2 , rather than simply escaping back 
into the gas phase. They argued that for gas and grain temperatures typical of molecular 
clouds, both S and f a should be of order unity. 

Hollenbach & Mckee (1979) later used this prescription to derive an H 2 formation rate 
for the local ISM that continues to be widely cited: 

A; 9 = 3xlO- 17 S^cm 3 s- 1 , (31) 
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where 



y = l + 0Ay/T 2 + T gr>2 + 0.2T 2 + 0.08T|, (32) 



and where T 2 and T grj2 are the gas and grain temperatures in units of 100 K. They argued 
that f a should be approximately constant and of order unity for grain temperatures below 
some critical value T cr , but that for T gr > T cr , it should fall off exponentially, with most of 
the hydrogen atoms evaporating from the grain surface before they have time to form H 2 . 
The value of T cr has proved hard to determine precisely, but is of the order of 100 K. 

Although this rate has been widely adopted in the literature, recent experiments have 
cast doubt on its validity at high temperatures, and suggest that the H 2 formation rate 
may be smaller than previously assumed (Pirronello et al. 1997a,b, 1999; Katz et al. 1999; 
Biham et al. 2001). However, since this conclusion is not entirely clear (Cazaux & Tielens 
2002) and their work is still ongoing, I have tentatively adopted the Hollenbach & Mckee 
rate below, with the proviso that the values of V cr that I derive may prove to be lower limits 
if the results of Pirronello et al. are borne out by future work. 



3. Comparing the different modes of formation 

We can combine equations 19 and 23 to write the total gas-phase H 2 formation rate as 

-Rfla.gas = -Rh 2 ,H- + -Rh 2 ,H+' ( 33 ) 

while the grain-catalyzed rate can be written as 

#H 2 ,dust = k 9 n tot n H \ ^—\ , (34) 

\L>MSN J 

where n tot is the total particle number density, V is the dust-to-gas ratio and Pmw is hs 
value in the local ISM. 

Combining these equations, we can easily solve for the dust-to-gas ratio at which i?H 2 ,gas 
and -RH 2 ,dust, which I denote as D cr : 

V cr = — — Pmw (35) 

In the common case that Rr 2 ,h- -Rr 2 > this equation reduces to 

T^cr = T^-p- / — - c — ^mwi (36) 

k 9 1 + x/x cr + / rad 
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Fig. 1. — V CT as a function of temperature for low ionization, low density gas (x = 10~ 4 , nn = 
1 cm" 3 ). The solid line is the value of V CT given by equation 35; the dashed and dotted lines 
are the contributions to this value of H and H^" respectively. 
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Fig. 2. — As figure 1, but for low ionization, high density gas (x = 10 4 , n H = 10 3 cm 3 ). 
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Fig. 3. — As figure 1, but for high ionization, low density gas (x = 10 2 , % = 1 cm 3 ). 




Fig. 4. — As figure 1, but for high ionization, high density gas (x = 10 2 , nu = 10 3 cm 3 ). 



where x cr and / r ad are given by equations 11 and 14 respectively. A similar equation can be 
written in the much less common case that R H2 H + -Rh 2 ,h-- 

In order to help illustrate the behaviour of these equations, I plot in figures 1 to 4 the 
value of V CT as a function of temperature for gas illuminated by the Mathis et al. (1983) 
radiation field in four different scenarios: low ionization, low density gas (x = 10~ 4 ,nn = 
1 cm -3 ; figure 1), low ionization, high density gas (x = 10 _4 ,nH = 10 3 cm~ 3 ; figure 2), high 
ionization, low density gas (x = 10~ 2 , nn = 1cm -3 ; figure 3) and high ionization, high density 
gas (x = lCT 2 ,nH = 10 3 cm -3 ; figure 4). In each case, I adopt a fixed grain temperature 
T gr = 20 K, although small changes in T gr have little effect on the results provided that it 
remains less than T cr . 

A striking feature of these plots is the strong temperature dependence of T> cv . At low 
temperatures, grain-catalyzed H 2 formation is relatively efficient and very little dust is needed 
before grain catalysis dominates. Above a few hundred K, however, the efficiency of grain 
catalysis decreases significantly, while the efficiency of gas-phase H 2 formation continues to 
grow. As a result, the required dust abundance rises sharply with increasing temperature. 



4. Astrophysical examples 

From the behaviour outlined in figures 1-4, it is clear that gas-phase H2 formation is 
at its most effective in warm, dense gas with a high fractional ionization. However, most 
of the molecular gas that we observe in our galaxy is in the form of molecular clouds with 
low temperatures (T ~ 20 K) and very low fractional ionizations (x ~ 10~ 7 ) and in these 
conditions grain catalyzed formation dominates by many orders of magnitude. 

A more promising place to look for gas-phase H2 formation is in the so-called warm 
neutral medium (WNM). In models of the multiphase ISM that assume thermal pressure 
equilibrium between phases (Field et al. 1969; McKee & Ostriker 1977; Wolfire et al. 1995), 
this is predicted to have a temperature of approximately 8000 K, high enough to collisionally 
dissociate H 2 . However, recent observations (Heiles & Troland 2002) and simulations that 
include the effects of turbulence (Gazol et al. 2001; Mac Low et al. 2001) suggest that much 
of this gas is actually at much lower temperatures; for instance, Heiles & Troland quote a 
temperature range of 500 < T < 5000 K. 

Taking representative values for the temperature and ionization of the WNM to be 
T = 2000 K and x = 10~ 2 (Heiles 2001), I find that 

V CI 4.1 0.047 
— — = 1 

X>MW 2.7 + / r ad,H- 11 + / ra d,H+ 



(37) 
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In other words, gas-phase formation would dominate if we could ignore the effects of the 
radiation field. In practice, this is not possible; at the densities characteristic of the WNM 
(n ~ 0.1 cm -3 ), we have 

/rad,H- = 3.9 XlO 3 (38) 

/rad,Ht = 1-0 XlO 2 , (39) 

and equation 37 becomes 

V CT = 5.2 x 10 -3 P M w- (40) 

These two examples demonstrate that gas-phase H 2 formation is unimportant in the 
bulk of the gas in the Milky Way: either the temperature and ionization are too low, as in 
molecular clouds, or the gas is too diffuse and H 2 formation is suppressed by the photodis- 
sociation of H~ and . 

However, there are a few counterexamples. For instance, gas-phase H 2 formation has 
long been known to play an important role in the chemistry of nova ejecta (Rawlings 1988) 
and protostellar outflows (Glassgold et al. 1989), where the gas initially has little or no 
dust (although more generally forms later) and where the high gas densities help mitigate 
the effects of photodissociation. Gas-phase formation is also predicted to dominate the 
molecular chemistry of freely-expanding supernova remnants such as SN1987A (Culhane 
& McCray 1995). Finally, Lepp & McCray (1983) suggest that gas-phase formation may 
dominate in X-ray dissociation regions (dense clouds illuminated by hard X-rays); modelling 
by Maloney et al. (1996) would appear to confirm this. 

These are somewhat unusual conditions, however, and in general dust abundances signif- 
icantly below the typical galactic value are required before gas-phase H 2 formation becomes 
competitive with grain-catalyzed formation. 

One place in which we might expect to find these low dust abundances is in the metal- 
poor gas within dwarf galaxies. Kamaya & Hirashita (2001) examine a well-studied example, 
the metal-poor dwarf IZwl8, and show that provided that its neutral ISM is clumpy (with 
clumps densities n > 100cm" 3 ) and moderately ionized (x ~ 10~ 3 ), then gas-phase formation 
will dominate. They also show that the formation of H 2 in this manner would not conflict 
with the upper limit on the H 2 column density of IZwl8 obtained by FUSE (Vidal-Madjar 
et al. 2000), due to the small filling factor of the clumps. 

Another place we might look for significant gas-phase H 2 formation is in damped Lyman- 
a (DLA) systems, many of which have low dust abundances (see, eg Lopez et al. 2002). An 
interesting example is the absorber at z — 3.025 in the spectrum of Q0347-3819 recently 
studied by Levshakov et al. (2002). This system has a temperature T ~ 400 K (as inferred 
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from the Doppler broadening of its many associated H 2 and metal absorption lines), and a 
fractional ionization i~2x 10~ 5 . If we assume that H~ photodetachment is negligible, we 
find that for this system 

D CI = 4.6 x 1(T 4 D MW . (41) 

Comparing this with a measured dust-to-gas ratio of T> ~ 0.05£>mw, we see that gas-phase 
H2 formation contributes no more than about 1% of the total H2 in this system. Including 
the effects of radiation merely strengthens this conclusion. 

There is no reason to suspect that this situation is particularly unusual; all damped 
Lyman-a systems by definition have large H I column densities, and consequently will have 
small fractional ionizations. We would therefore expect grain-catalysis to dominate in these 
systems. 

Finally, gas-phase H 2 formation has long been known to play an important role in the 
early stages of galaxy formation. In primordial gas, this is obvious: there is no dust, so any 
H2 that forms must form in the gas phase. A more interesting problem is determining the 
value of T> cv for these systems; in other words, at what point does grain- catalyzed formation 
overtake gas-phase formation? 

For the purposes of this discussion, I adopt the example of an H 2 -cooled protogalaxy 
with temperature T = 1000 K and fractional ionization x = 2 x 10~ 4 (Tegmark et al. 1997). 
These values are appropriate for the first generation of star-forming protogalaxies, and while 
they may be underestimates for later generations, my analysis can easily be rescaled for 
higher values. For this example protogalaxy, I find that 

V CT 2.1 x 10~ 2 1.3 x 10~ 4 
— — = 1 

^MW 1 + /rad,H- 1 + / ra d,H+ 

Thus, if the radiation field is unimportant, V CT ~ 0.02 T>mw, comparable to the values seen 
in some metal-poor dwarf galaxies at the present-day (Lisenfeld & Ferrara 1998). 

How strong is the radiation field within a protogalaxy? There are potentially three 
main contributors to this field: the protogalaxy's own stellar population, emission from 
neighbouring galaxies and radiation from the cosmological background produced by distant 
sources. It is simplest to consider these separately. 

Much of the optical and near-infrared radiation responsible for destroying H~ and is 
produced by long-lived stars, and so the contribution of the protogalaxy's stellar population 
depends as much on its star formation history as on its current star formation rate. This 
makes it very difficult to parameterize its effects in the general case; it is much easier to 
examine a simple example that will hopefully be broadly representative. 



(42) 
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For the purposes of this example, I assume: 

1. That the protogalaxy underwent an instantaneous (or near-instantaneous) starburst 
10 8 yr ago, following which it has formed no more stars. 

2. That the stars which did form are located in the centre of the protogalaxy, within a 
small enough region that I can approximate their emission as coming from a point 



3. That the luminosity and spectral energy distribution of this stellar cluster are well 
described by the Z = 0.05Z Q model of Leitherer et al. (1999). 

4. That the protogalaxy itself is well-described by a truncated isothermal sphere density 
profile (Iliev & Shapiro 2001). 

All of these assumptions are debatable, but they do provide us a basis on which to estimate 
the effects of the stellar radiation field. Moreover, these assumptions are somewhat conser- 
vative, and tend to minimize the effectiveness of the stellar radiation. For instance, if we 
reduce the time since the starburst from 10 8 yr to 10 7 yr, then the photo detachment rate 
increases by a factor of fifty. Similarly, if we assume continuous star-formation rather than 
an instantaneous starburst, then we obtain a similar (or slightly larger) photodetachment 
rate once the total mass of stars formed has reached a comparable level. 

We could also criticize the adoption of the Leitherer et al. model, on the basis that it 
assumes a standard Salpeter IMF, while there is considerable evidence that the the primordial 
IMF is biased towards high masses (Larson 1998). However, this again means that we will 
underestimate the photodetachment rate (although we will significantly overestimate the 
lifetime of the stellar population). 

Returning to my example, the first two assumptions allow me to write the H photode- 
tachment rate at a distance R from the stars as 



where L„, the stellar luminosity per unit frequency, is given by the Leitherer et al. model. 
Using this value, I obtain 



source. 




(43) 



k 6 = 4.9 x 1Q- S F(M^R) s 



-i 



(44) 



where 




(45) 
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and where M* is the mass of stars formed in the starburst. Similarly, we can write the 
photodissociation rate as 

k 8 = 9.4 x 10~ 10 F(M*, R) s~\ (46) 
and from these rates calculate / ra d,H- and / radjH + 

/rad,H- = AOn^F(M^R) (47) 
/rad,H+ = 1.5rgF(M.,R). (48) 

By comparing these values and equation 42, we can see that formation via contributes 
at most about 10% of the H 2 produced in the gas phase, with the rest coming from H~. 

To evaluate these numbers, I use the fact that for a truncated isothermal sphere, 

R \ „ „ _ „ q 



n H _— ~6.9xl0- 3 -^ (49) 

Vikpc/ n m 

in regions outside of the core. The final unknown, M*, can be written as 

M* = 1.5 x 10 7 e*(l + z)- 3 / 2 ( ] M , (50) 

where e* is the star formation efficiency of the protogalaxy, z is its redshift of formation and 
h is the Hubble constant in units of 100 km s _1 Mpc -1 . 

For a protogalaxy that formed in a standard ACDM cosmology (Q m = 0.3, ilb = 0.04, 
h = 0.7) at a redshift z = 10, and that formed stars with an efficiency a* = 0.01, we find 
that 

V CI ~ 3.8 x 10~ 4 P MW , (51) 

Thus, in this particular example, radiation from the existing stellar population reduces D cr 
by almost two orders of magnitude. 

In view of the uncertainties involved in producing this estimate, it would be unwise 
to over- generalize. However, since my assumptions verge on the conservative side, it seems 
likely that in realistic protogalactic models we would see similar effects, and that gas-phase 
H 2 formation will rapidly be overtaken by grain- catalyzed formation. 

What about protogalaxies that have yet to form stars? In this case, there is no significant 
local contribution to the radiation field, which instead is produced by neighbouring sources 
and/or the cosmological background. 

For neighbouring sources, we can reuse the above formalism, as long as we set R to 
the distance to the extragalactic source. However, this is typically an order of magnitude 
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or more greater than the size of a protogalaxy, implying that the effect of the radiation will 
be at least two orders of magnitude smaller than the effects discussed above. Consequently, 
radiation from protogalaxies of the size discussed here will have little or no effect on gas- 
phase H 2 formation within their neighbours, unless their emitted flux is substantially larger 
than has been assumed here. 

For the background, we again face the problem that any conclusions that we can draw 
are strictly limited by our poor knowledge of the star formation history, this time on a 
cosmological rather than protogalactic scale. The best that we can do is to determine how 
strong the background needs to be before it has a significant effect. Modeling the background 
below the Lyman limit as a power-law, 

j»=j2i(^y, (52) 

where J21 = 10~ 21 erg s" 1 cm~ 2 Hz -1 sr^ 1 and where hv = 13.6 eV, I find that for a = 1, 

k 6 = 8.1 x l(T 10 J 2 i s _1 , (53) 

and hence 

/rad,H- =0.6 J 21 n R \ (54) 

The significance of the background varies with Uh and hence with position within the 
protogalaxy. For my example protogalaxy, formed at a redshift z — 10, gas near the trun- 
cation radius has a density n H ~ 7 x 1CT 3 cm~ 3 and thus is affected for J 2 i ^ 0.01; on the 
other hand, gas in the central core has % ~ 1 cm~ 3 and is only affected for J21 > 1. 



5. Conclusions 

The simplicity of the basic chemistry involved in gas-phase H2 formation means that it 
is easy to construct a fairly accurate expression for the formation rate in terms of only a few 
parameters: the temperature, density and fractional ionization of the gas, plus the strength 
of the radiation field near the H~ and H_2 photodissociation thresholds. Expressions for the 
rate of H 2 formation via the H and ions are given by equations 19 & 23 respectively, 
and the total formation rate is simply the sum of these two values. 

Using these expressions, together with an analytical expression for the grain- catalyzed 
H 2 formation rate, one can solve for V cr , the dust-to-gas ratio required for grain- catalyzed 
H_2 formation to overtake gas-phase formation. The results demonstrate that, in principle, 
gas-phase H 2 formation could be comparable to grain-catalyzed formation in galactic gas, 
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particularly at high temperatures where the latter is inefficient. In practice, however, it is 
usually significantly slower, either because of a shortage of free electrons and protons (which 
reduces the formation rate of the intermediate ions) or because the incident radiation field 
destroys the ions before they have a chance to form H 2 . 

Finally, in order to demonstrate the simplicity and potential utility of this approach, 
I have applied it to a number of astrophysical examples. Not surprisingly, I find that in 
most cases T> CT is significantly less than the mean galactic value, often by several orders of 
magnitude. Nevertheless, there are counterexamples, such as X-ray photodissociation regions 
(Lepp & McCray 1983; Maloney et al. 1996), nova ejecta (Rawlings 1988), or the high-redshift 
protogalaxies analyzed in detail here. A common thread linking many of these exceptions 
seems to be the fact that they have dust-to-gas ratios (but not necessarily metallicities) 
significantly lower than the mean galactic value. 

I would like to acknowledge useful comments on an earlier draft of this paper from 
Mordecai Mac-Low and Michael D. Smith. I would also like to thank the anonymous referee 
for a timely and useful report, and for bringing the work of Lepp & McCray to my attention. 
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Table 1. Reaction rates 



Reaction Rate (cm s ) Reference 

1. H + c^H-+7 h = 1.4 x 10~ 18 T a928 exp (-^o) de Jong (1972) 



n TT , n , i 1.5 x 10- 9 T< 300K , , ,^ nM . 

2. H + H - H 2 + e k 2 = { 4 q x 10 _ 9T _ . 17 y 3Q0K Lannay et a/. (1991) 



^ t t TT i TT ^ , dexf— 19.38 — 1.523 log T , . , ,„„_„ x 

3. H + H+ H+ + 7 * 9 = + l imog(Tf _ .4 91og(r)3] Ramaker & Peek (1976) 

4. H+ + H -> H 2 + H+ A; 4 = 6.4 x lO" 10 Karpas e* a/. (1979) 

r TT tt+ nTT , 5.7 x i(r 6 T- a5 + 6.3 x 1(T 8 , 

5. H- + H+^2H k 5 = 11mn _ A Moseley e* a/. (1970) 

5 -9.2 x 10- n T a5 + 4.4 x 10~ 13 T J v ; 

6. H + 7 — > H + e See text 

n. TT + ott , f 1.0 xlO" 8 T<617K _ . . 

7. H+ + e^2H A; 7 = ^ „ _ 7fi _ _ Schneider e* a/. (1994) 

2 7 \ 1.32 x Kr 6 T-°- 76 T > 617K v ; 

8. H+ + 7 -> H + H+ See text 

9. H + H + grain — * H 2 + grain See text 



